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Abstract 



We present a phenomenological discussion of spin-independent, general- 
ized dipole polarizabilities of hadrons entering the virtual Compton scattering 
process j*h — > jh. We introduce a new method of obtaining a tensor basis 
with appropriate Lorentz-invariant amplitudes which are free from kinemat- 
<D \ ical singularities and constraints. The result is summarized in terms of a 

compact effective Lagrangian. We then motivate a gauge-invariant separation 
into a generalized Born term containing ground-state properties only, and a 
^ ■ residual contribution describing the model-dependent internal structure. The 

generalized dipole polarizabilities are defined in terms of Lorentz-invariant 
residual amplitudes. Particular emphasis is laid on a physical interpretation 
of these quantities as characterizing the spatial distributions of the induced 
electric polarization and magnetization of hadrons. It is argued that three 
dipole polarizabilities, namely the longitudinal electric a^^q 2 ), the transverse 
electric axiq 2 ), and the magnetic (3(q 2 ) ones are required in order to fully re- 
construct local polarizations induced by soft external fields in a hadron. One 
of these polarizabilities, ar(q 2 ), describes an effect of higher order in the soft 
final-photon momentum q' . We argue that the associated spatial distributions 
obtained via the Fourier transforms in the Breit frame are meaningful even 
for such a light particle as the pion. The spatial distributions are determined 
at large distances r ~ l/m^ for pions, kaons, and octet baryons by use of 
ChPT. 
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I. INTRODUCTION 



Recently, Compton scattering with virtual photons (virtual Compton scattering, or VCS) 
has attracted considerable experimental and theoretical interest (for an overview see, e.g., 
]T|-[3|] and references therein). Both the near-threshold and high-momentum transfer 
regimes f^|7| of VCS turn out to be very informative for studying the structure of hadrons. 

In the below-threshold region a set of structure functions, the so-called generalized polar- 
izabilities (GPs), has been introduced in Ref. @ in order to parametrize structure-dependent 
effects in the VCS amplitude to leading order in the final-photon momentum q'. These GPs 
contain model-dependent information beyond the low-energy theorem (LET) of real |||| 
and virtual Compton scattering |IU |. The first experimental results for generalized polariz- 



abilities of the proton have recently been obtained at the Mainz Microtron (MAMI) for a 
four-momentum squared of Q 2 = 0.33 GeV 2 [IT]. Additional experiments aiming at proton 



polarizabilities are presently carried out at Jefferson Lab [12| and MIT-Bates fL3|. A sen- 



sitivity study of inelastic high-energy pion-electron scattering on the pion VCS amplitude 



has been performed as part of the Fermilab E781 SELEX experiment [[L4],|15|]. 

Theoretical predictions of the nucleon GPs have been obtained in different frameworks 
such as the constituent quark model |pj|16|-^8[|, the linear cx-model [jlJJ, chiral perturbation 



theory (ChPT) pIETI], the Skyrme model |]22||, models with effective photon-pion- nucleon 



Lagrangians |23|,|24[], and dispersion theories p5|j26[1 . Pion and kaon GPs have been discussed 
at the one-loop level of ChPT |Z7||Z§, 

In the present paper we are not so much concerned with specific results for GPs as 
obtained from various models of hadrons. Our investigations are rather directed to an 
analysis of general properties of GPs as well as their physical interpretation which, in our 
opinion, was only insufficiently covered previously. We mainly consider kinematical aspects, 
relegating dynamical aspects of properties of the generalized polarizabilities such as sum 
rules, dispersion relations, etc. to a future publication. In order to avoid complications 
related with the spin of the target we only discuss the simplest case of a (pseudo) scalar 
particle like the pion. Our considerations may as well be applied to the spin-averaged part 
of the VCS amplitude for a target of arbitrary spin. For illustrative purposes we make use 
of results obtained in the framework of ChPT. 

Our work is organized as follows. In Sec. II, we introduce a new method of obtaining 
invariant amplitudes of electromagnetic reactions free from kinematical singularities and 
constraints, using (double) virtual Compton scattering as an example. Section III contains 
a motivation for choosing a specific (though standard) form of the Born amplitude which 
is used for separating convectional from internal contributions in the scattering amplitude. 
In Sec. IV, we show how generalized dipole polarizabilities can be defined in a Lorentz- 
invariant manner from the invariant amplitudes without introducing inconvenient singular 
kinematical factors. In particular, we show that it is natural to introduce one more dipole 
polar izability, namely, the transverse electric one, which is needed in order to recover the 
electric polarization of the target induced by a soft external electric field. This quantity does 
not contribute in the 0(q') limit of previous analyses. We provide another explanation why 
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to this order only two of the three spin-independent polarizabilities are accessible in photon 
electroproduction experiments. We give a physical interpretation of the generalized dipole 
polarizabilities in terms of densities of the induced electric and magnetic polarizations. In 
Sec. V we argue that the associated space distributions are meaningful even for a light 
particle like the pion. We determine such space distributions at large distances r ~ l/m n 
derived from various form factors using the predictions of ChPT. In passing, we also include 
SU(3) / extensions of previous results obtained within two- flavor HBChPT. Analytical results 
of spatial distributions obtained via dispersion relations are relegated to the Appendix. 



II. TENSOR BASIS AND INVARIANT AMPLITUDES OF VCS 

In our discussion of virtual Compton scattering off a spinless hadron, we will often refer 
to this hadron as a "pion," although our general considerations also hold true for other 
spinless particles, nuclei, and even atoms, as well as for spin-averaged amplitudes. We 
start our analysis with an investigation of the general kinematical structure of the pion VCS 
amplitude Tycs for the case of two virtual photons, assuming that both initial and final pions 
are on shell. Our aim is to construct a tensor basis and a related set of Lorentz-invariant 
amplitudes B iy that provides a complete parametrization of Tycs- We require all Bi to be 
free from kinematical singularities and constraints, because this simplifies the classification 
of low-energy characteristics of the pion, and also provides technical advantages, for instance 
when discussing dispersion relations. 

The problem of finding a set of amplitudes for VCS was already addressed and solved by 
Tarrach for both spin-0 and spin-1/2 targets by using a projection technique originally 
proposed by Bardeen and Tung |30[] . Although, in principle, we could directly use the results 
of Ref. P?l, we prefer to construct the tensor basis again, first, in order to demonstrate a very 
simple and powerful alternative method which avoids projections and, second, in order to 
rearrange the VCS tensor in a manifestly gauge-invariant form and to split it into structures 
contributing to real Compton scattering, VCS with one photon virtual and, finally, VCS 
with both photons virtual. 

To begin with, we define the amplitude Tycs of virtual Compton scattering, 



where T^ v is the Compton tensor given in terms of the covariant (Wick) Ty^-product of the 
electromagnetic currents, 




(2.1) 



as 



Tvcs = e^T^, 



(2.2) 




(2.3) 
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We normalize all single-particle states as 

(p'\p) = 2p (2n) 3 5 3 (p-p'), (2.4) 

so that in case of real photons the S'-matrix reads 

S ft = t(2n) 4 5 4 (p + q - p >-q>) T ycs . (2.5) 

Due to four-momentum conservation, p+q = p' + q', the tensor depends on three linearly 
independent vectors P, Q, and R: 

p = \(p + p), Q = \(q + q), R = \{p -p) = \{q-q')- (2.6) 

Since we only consider the case of initial and final pions on the mass shell, the vectors P, 
Q, and R are constrained by 

P 2 = M 2 — R 2 , P-R = 0, (2.7) 

where M is the pion mass. Hence, we can choose four independent kinematical invariants 
which, for the moment, we take to be Q 2 , R 2 , P ■ Q, and Q ■ R. 

The discrete symmetries as well as gauge invariance impose restrictions on the general 
form of the Compton tensor. For charged pions the combination of pion crossing with 
charge-conjugation symmetry results in[] 

T^(P, Q, R) = T^i-P, Q, R), (2.8) 

whereas photon crossing yields 

7V(P, Q, R) = T UIM (P, -Q, R). (2.9) 



Gauge invariance requires 



= (Q + R)^ = 0, 
q' v T^ = (Q - R) V T^ = 0, (2.10) 



where the first and second equation are not independent once photon-crossing symmetry is 
imposed. Finding a solution to Eqs. ( [2.10| ) without introducing kinematical singularities or 
constraints is the main challenge in constructing appropriate amplitudes. 



1 For the neutral pion (but not for the neutral kaon), which is its own antiparticle, only pion 
crossing is required to obtain Eq. (|2.8|), 
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Because of parity conservation the Compton tensor transforms as a proper second-rank 
Lorentz tensor. The most general T^ u satisfying the crossing-symmetry conditions of Eqs. 



[q) and ( |2.9| ) can be written as a linear combination 

10 



i=l 



of ten basis tensors r* which include the metric tensor g^ v and 9 bi-linear products of P, 
K, and Q:Q 



r 1 
1 /if 


fi'/if > 




r 2 
/if 


= P P 




r 3 
1 /tf 


QiiQvi 




r 4 
1 /tf 






r 5 
/tf 


= (P tl Q u + P„Q l *)(P- 


Q), 


r 6 
/if 


= (P^Q„-P.Qv)(P 


Q){Q 


r 7 
1 /if 


= (P^R V + P V R^)(P ■ 


Q)(Q 


r 8 
/if 


= (P^R U - P U RJ(P ■ 


Q), 


r 9 
/if 


= {Q^R U + QuRfi){Q 


■R), 


r 10 
1 /if 


= (Q/iRv — QvR/jt)- 





(2.12) 

Appropriate factors of P ■ Q have been introduced such that all r' are even functions of P, 
as required by Eq. Q2.8p . Correspondingly, factors of Q ■ R provide for photon-crossing-even 
basis tensors [see Eq. ( p.9|) ]. With such a choice of t*, all the functions Ai depend on the 
crossing-even variables Q 2 , R 2 , (P ■ Q) 2 , and (Q ■ R) . 



At this point it might be worthwhile to explain why introducing factors of P ■ Q and Q ■ R 
into Eqs. ( 2.12 ) is harmless to the analytical properties of the functions Ai in Eq. ( |2.11| ). For 



that purpose, let us omit in Eqs. ( |2.12|) all factors of P ■ Q or Q ■ R, and denote the resulting 

o 

basis tensors by r* with Ai the corresponding functions of the expansion of T^ v in terms of 
the "reduced" tensors f* Let us further assume that the tensor results from some set of 
Feynman diagrams consistent with all symmetries. Consider an arbitrary Feynman diagram 
denoted by G. Obviously, its contribution to T^ u can be expressed in terms of or bi-linear 
products of four-momenta (these are just the r % ) multiplied by scalar coefficients. Any such 
coefficient will at most have dynamical singularities related with propagators of intermediate 

o 

particles but no kinematical singularities. Stated differently, all individual contributions Af 

o 

to the functions Ai are free from kinematical singularities. In general, a single contribution 
is not separately crossing symmetric. Eventually, crossing symmetry of is obtained after 



2 Due to parity, structures containing a single fully antisymmetric tensor are excluded. 
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adding one or several crossed partners, G c , of the diagram G. Let us consider, for instance, 
the tensor structure = P^Q U + P V Q^ which is odd under P — > —P. We can then write 

A G+G ^x)=A G (x)-A G (-x), (2.13) 

where we introduced x = P ■ Q and, for brevity, omitted P-independent arguments like 
Q 2 . The second term in Eq. ( |2.13| ) represents the contribution of the crossed diagram G c 

which makes the amplitude A G+Gc (x) odd. Note that both contributions are non-singular 

as x -+ 0. From Eq. (12~T3] ) we conclude that A G+Gc (x) = xA G+Gc (x), where A G+Gc (x) is 
an even function of x which has no pole at x = and which therefore has no kinematical 
singularities at all. In other words, the sum of all Feynman diagrams contains the term 
r^Ac,, the tensor r£, carrying a factor of P ■ Q and the function A§ having no kinematical 
singularities. 

The same consideration immediately applies to the tensors t^ v and t^ u . In the case of 
or tI v , we have to apply the procedure given by Eq. ( |2.13| ) twice — first for showing 

o 

that the corresponding function Ai contains the factor of P ■ Q, and second for showing that 
it contains the factor of Q ■ R as well. Eventually, we conclude that all the coefficients Aj in 
Eq. ( |2.11 ) are free from kinematical singularities. 



By means of the factors P ■ Q and Q ■ R in the tensor basis Eq. (|2.12|) we can solve the 



constraints of Eqs. (|2.10| ) without introducing singular coefficients and without using the 
projectors suggested in Ref. Indeed, inserting Eq. ( [2.11|) for T^ v into Eqs. ( j2.10|) and 



collecting coefficients of the independent four-momenta P, Q, and R, we obtain a set of six 
linear equations in the functions A it of which only five are independent: 

A 2 + Q 2 A 5 - (Q ■ R) 2 A 6 + (Q ■ R) 2 A 7 - R 2 A 8 = 0, 

A 5 - Q 2 A 6 + R 2 A 7 - A 8 = 0, 

A, + Q 2 A 3 + {P ■ Q) 2 A 5 + (Q ■ R) 2 A 9 - R 2 A W = 0, 

A 3 + (P • Q) 2 A 6 + R 2 A 9 - A w = 0, 

A, + R 2 A 4 + (P • Q) 2 A 8 + (Q ■ R) 2 A 9 + Q 2 A 10 = 0, 

A 4 + (P ■ Q?A 7 + Q 2 A 9 + A 10 = 0. (2.14) 

One can now express Ai, i — 1, ■ ■ ■ , 5, in terms of the remaining functions Ai, i — 6, • • • , 10, 
without singular coefficients: 



l 10; 



A 1 = R 2 (P ■ QfA 7 - (P ■ Q) 2 A S + [Q 2 R 2 - (Q ■ R) 2 ]A 9 + (R 2 - Q 2 )A, 
A 2 = l(Q ■ R) 2 - (Q 2 ) 2 }A 6 + [Q 2 R 2 - (Q ■ R) 2 }A 7 + (P 2 - Q 2 )A 8 , 
A 3 = -(P-Q) 2 A 6 -R 2 A 9 + A 10 , 
A 4 = -(P-Q) 2 A 7 -Q 2 A 9 -A 10 , 

A 5 = Q 2 A 6 -R 2 A 7 + A 8 . (2.15) 



Using Eqs. (|2.15| ), we can finally rewrite Eq. ( [2.11| ) as 
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T,u = E T; U A(Q 2 , R 2 , (P ■ Q)\ (Q • R) 2 ), (2.16) 



i=6 



where 



7* = [(Q • J2)^ - {Q'Y]P,P U - (P ■ QYQ^Q, 



+ Q 2 (P ■ Q)(P^Q U + P U Q,) + (P ■ Q)(Q • R)(P»Q U - P V QJ, 
1%, = R\P ■ Q) V + IQ 2r2 ~ (Q ■ Rf\P^ ~ (P ■ QfR^Ru 
-R 2 (P- Q)(P,Q V + P U Q,) + (P • Q)(Q ■ R)(P,R U + F„ig, 
= -(P ■ Q) V + [F 2 - Q 2 \P,P V 

+ (p • q)(p„q v + P U Q,) + (P • Q){P,R V - P„RJ, 

T^u = [Q 2 R 2 — (Q ■ R) 2 ]9tiv — R 2 Qfj,Qu — Q 2 R^R U + (Q ■ R)(Q^Ru + QuRfi), 

TH = {R 2 - Q 2 )g^ + Q^Q U - R^R, + (Q^ - Q v R p ) (2.17) 

are five basis tensors which explicitly satisfy the crossing-symmetry and gauge-invariance 
conditions of Eqs. (|2.8|)-( p7L0|) . Respectively, the five functions Ai, i = 6, • ■ ■ , 10, can be 
considered as invariant amplitudes of virtual Compton scattering which are free from kine- 
matical singularities and constraints. 

In passing we note that the same method of constructing a basis and invariant amplitudes 
free from kinematical singularities and constraints also works for the VCS amplitude of a 
spin-1/2 target such as the nucleon. 

The tensors T®„ and have exactly the same form as for real Compton scattering and 
can easily be identified with the more common notation 

1%, = ~(P -Q)(P- Q')9^ - (? ■ 4)P»Pv + (P ■ q')P,q, + (P ■ q)P^, 

TH = %qu-{q-q)g,u- (2.18) 



However, also the remaining tensors T l v , i = 6, 7, 9, and thus the corresponding functions 
Ai contribute to RCS. This, clearly, is a drawback of the basis of Eqs. ( |2.17| ) and it would 
be convenient to have, instead, a basis such that exclusively the two tensors T^ u and 
contribute to RCS, another one appears for the case of one virtual photon and, finally, 
the two remaining structures also contribute to 7*7r — ► 7*7r. To that end let us introduce 
gauge-invariant combinations of photon polarizations and momenta, 

F ia , = -i(q l *e v -q v e ll ), F%, = itftf - <£ej). (2.19) 

These second-rank tensors represent the Fourier components of the electromagnetic field- 
strength tensor J-^x) = d^A^x) — dyA^x) associated with plane-wave initial and final 
photons described by vector potentials A^(x) = e M exp(— iq ■ x) and A' (x) = e'* exp(iq' ■ x), 
respectively. In terms of F^ v and F' it turns out to be rather straightforward to identify 
structures contributing for real or virtual photons. 
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For example, the A w contribution to the VCS amplitude reads 

e»e'*%°A 10 = -If^F'^Aw (2.20) 

which can be interpreted as the matrix element of the effective LagrangianP] 

£ = -^10 V^V- (2-21) 

Here A\o represents a differential operator in terms of (covariant) derivatives acting on 
both pion and electromagnetic fields with Fourier components given by the function Aiq. 
Similarly, the contribution of the amplitude A$ can be written as 

e^T^As = -(P,Fn(P p F' pi/ )As (2.22) 

which results from the effective interaction Lagrangian 

C = ~A^ av T Pv P a P p ^(j), (2.23) 

where the action of PaP^cj) is defined by 
P^<P= l -^D p <p- l -(D^^ 
P a P^(j> = --tfD a Dp<f> + j{D a (j)^Dp(t> + ^{D p( j^D a (j> - ^{D a D P (j^(j>, (2.24) 

and D a cj) = d a (f) + ieZA a cj) with e > 0, e 2 /47r « 1/137, and Ze denoting the charge of the 
particle. 

The tensor structures T l v , i = 6, 7, 9, involve higher powers of photon momenta and thus 
are related with derivatives of the electromagnetic fields. Introducing the four-vectors 

- iq^F^ = -q 2 e v + (q ■ e)q v , 

iq >»F^ = -q' 2 e'; + (q'.e'*)C (2-25) 

which vanish for real photons, we obtain two more tensors, namely T® and T® — T^, by 
using the identities 



3 Here, we assume that "pion" and "antipion" are different particles, such as 7r + and tt or K° and 
K°. For the case of a charged pion, the field (p = 7r+ is given in terms of the Hermitian, Cartesian 
isospin components 4>i as <fi = (4>i — i</>2)/V2 (4 1 = vr~ = {4>i + i<p2)/V^) and destroys a ir + (vr~). 
In the case of the neutral pion, we have to take <p = <p^ = 4>3 and replace the factor of 1/4 in Eq. 
( 2.21| ) by 1/8. These trivial changes also apply to other Lagrangians written below. 
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iMtFnti'K") = ^* V \.K» + ( r2 - Q 2 K°i ( 2 - 26 ) 

and 

(P„ qtl Fn(P a q' P F' pa ) = eV"[-27j, + 2l£, + (i? 2 + Q 2 )T«, - (P ■ Q?T$. (2.27) 



The remaining linear combination T pu + Tj^ is contained in a product similar to Eq. ( p.27|) 
however with photon momenta interchanged in one factor: 

(P v F^q')(PW p F') + {P v F> q »){P a q p Fn 



pa J ' \ fjtu - 

eV^-ZZj, - 2Tl - 2R 2 Tp- (2-28) 



With the above identities the most general VCS amplitude can be written in the following 
manifestly gauge-invariant form: 

Tvcs = \f^F^B x + (P,Fn(P p F' pu )B 2 

+ [(P»q'»F pu )(P°q'»F' pa ) + (P»q»F^)(P°q<>F> a )]B 3 

+ (q,FnW P F' pu )B, + (P»q»F p „)(P°q'CF' pa )B 5 . (2.29) 

Here all the invariant amplitudes Bi are free from kinematical singularities and constraints, 
because the transformation from the basis T l of Eq. (|2.17 ) to the basis of Eq. (|2.29|) is 



non-singular. This can also be easily seen from the following relations between the two sets 
of amplitudes A4 and Bf. 

Bi = -An + - A T ) + ^=^4>, 
B 2 = -A 8 + —{A, + A 7 ) j^~(A e - A 7 ), 

B 3 = ~(A (i + A 7 ), B, = ^A 9 , B 5 = j(A 7 -A 6 ). (2.30) 

The determinant of the transformation expressing Bi (i = l,---,5) in terms of A$ (i = 
6, ■ • • , 10) is 1/32 7^ 0, independently of the values of the kinematical variables. 



In view of the rather compact and transparent structure of Eq. ( |2.29|) , we will use in the 



following the parametrization of Tycs given by the amplitudes B^. As seen from Eqs. ( |2.25|) 
and ( |2.29|) , only the amplitudes Bi and B 2 are needed to describe real Compton scattering, 
because then q^F^ = q'^F' = 0. When only one photon is virtual, one more amplitude 
(B 3 ) contributes. All five amplitudes Bi enter, when both photons are virtual. 



Equation fl2.29p can be interpreted as the matrix element of the effective interaction 

1 

4 

1 



C = -[B 1 {F pu ) 2 + 2S 4 (^J>) 



+ -\B 2 T^Tp p + B^d^WT^) - 2B 3 T^[d p d v T Pv )\P a P%^, (2.31) 



9 



where Bi are differential operators acting on all the fields and determined by their Fourier 
components Bi. 

Of course, after substituting q — > q x , q' — > —q 2 , p — > — pi, and p' — > P2, Eq. ( |2.29| ) 
also describes the general kinematical structure of the amplitude of the crossed reaction 
7(91)7(92) - * 7r (Pi) 7r (P2) for on-shell pions. Exactly the same considerations apply to any 
other spin-0 hadron assuming the same symmetry principles (Lorentz and gauge invariance, 
P, T, and C conservation) and are also applicable to a properly spin-averaged VCS amplitude 
for hadrons with finite spin. 

As mentioned before, the functions Bi depend on the four invariants Q 2 , R 2 , (Q ■ R) 2 
and (P ■ Q) 2 . As an alternative, the following combinations of the first 3 quantities can 
be used as independent arguments of Bf q 2 + q' 2 = 2(R 2 + Q 2 ), q ■ q' = R 2 — Q 2 , and 
q 2 q' 2 = (# 2 + Q 2 ) 2 - 4(Q ■ R) 2 . Thus, we may write 

B i = B i (v 2 ,q-q' : q 2 + q' 2 ,q 2 q' 2 ), (2.32) 

where v is defined as 

Mv = P-Q = P-q = P-q'. (2.33) 

Besides being manifestly crossing symmetric, this form of Bi has the advantage of having a 
simple limit if one or both photons become real. 

Finally, the Mandelstam invariants of the VCS reaction read 

8 = (p + q) 2 = M 2 + 2Mu + q-q', 
U=(p- q'f = M 2 — 2Mv + q-q, 

t=(q~ q'f = q 2 + q 2 ~ 2q ■ q. (2.34) 



III. THE BORN TERMS AND GAUGE INVARIANCE 

In order to describe the internal structure of the pion in terms of its generalized po- 
larizabilities as tested in virtual Compton scattering, we first have to isolate a convection 
contribution which originates in two successive interactions of the photons with the elec- 
tromagnetic current of the pion, resulting in singularities at zero photon momenta. For a 
point-like (pseudo-) scalar particle of electric charge eZ, the interaction with an external 
electromagnetic field can be described in terms of the Lagrangian[] 



4 Since we do not treat the electromagnetic field as a dynamical variable, it will not be included 
in the list of arguments of the Lagrangian. 
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where the covariant derivative 

D fi<f> = {d„ + ieZA„)(f> (3.2) 

ensures the invariance of the Lagrangian under the canonical gauge transformation 

Ap(x) I— > Afj,(x) + dfj,A(x), (f>(x) i— > exp[— ieZA(x)]<f>(x). (3.3) 

However, such a description is not sufficient for an extended particle and we have to modify 
the above effective Lagrangian. 

To that end, let us first consider a classical system of n constituents with charges e a = eZ a 
and masses m a which is exposed to a static external potential A (r). The electrostatic energy 
of such a system is given by 

n n i 

W = J2 e a A (r a ) = e a [A {R) + p ai ViA (R) + -p a4 p aj V i V j A ( J R) + •••], (3.4) 

a=l a=l ^ 

where R is the center of mass of the charge distribution and p a = r a — R are relative 
coordinates. In the continuum limit, the expression for a spherically symmetric distribution 
reads 

W = eZA (R) + ^(Zr 2 E )V 2 A (R) + ■■■ = eF(V 2 )A (R), (3.5) 
o 



where J2a e a — » / p{r) dr = Ze is the total charge, J2a e aP 2 l — ► J r 2 p{r)dr = e{Zr 2 E ) is 
the electric mean square radius, and F(—q 2 ) = Z — ^{Zr 2 E )q 2 + • • • is the electric form 
factor. Of course, the response of the extended system to the external field is determined 
by the potential and its derivatives together with the corresponding moments of the charge 
distribution. 

A relativistic generalization to the case of an extended pion suggests the following sub- 
stitution for the vector potential in Eq. ( |3.2| ) 

ZA^x) -> F(-d 2 )A^x) (3.6) 

which leads to the effective LagrangianP] 

C ee (d^, 0) = [d, + ieF(-d 2 )A,]<P[d» - ieF(-d 2 )A^ - M 2 0, (3.7) 



Since we want to apply the Lagrangian for the case of space-like virtual photons, we assume the 
form factor to be real. 
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and to the first-order electromagnetic vertex 

r») = (p' +p)»F(q 2 ), q = p'-p, F(0) = Z. 



(3.8) 



Note that Eq. (|3.7|) is no longer invariant under the canonical gauge transformation of 
Eq. (|3.3|) and, correspondingly, the electromagnetic vertex of Eq. (|3.8j) does not satisfy the 
Ward-Fradkin-Takahashi identity |32|-g4]] 

q ^(p',p)^Z[A~\p')-A-\p)}, (3.9) 

where A(p) = l/(p 2 — M 2 ) is the free propagator of C (d^(p, 0). In fact, a different gauge 
transformation 

A^x) i-> A^x) + <9 M A(x), (f>{x) i-> exp[-ieF(-<9 2 )A(x)]0(x), (3.10) 

defines a local realization of the symmetry group U(l), i.e., under two successive transfor- 
mations described by smoothly varying functions Ai and A2, the fields transform as 

A^ 1— > A^ + <9 M Ai i-> (A M + ^Ai) + «9 M A 2 = A^ + 9 M (Ai + A 2 ), 
<p h+ exp[-ieF(-cf i )A 1 ]<f> h-> exp[-ieF(-d 2 )A 2 ] exp[-ieF(-<9 2 )A 1 ]0 

= exp[-2eF(-9 2 )(A 1 + A 2 )]0. (3.11) 

Accordingly, we define a noncanonical covariant derivative as 

f>A= [d» + ieF(-d 2 )A^ (3.12) 

such that D^cj) transforms in the same way as under Eq. ( |3.10| ) and the effective Lagrangian 
C eS (d^, 0) = C (D^, 0) = D^iD^Y - M 2 00t (3.13) 

remains invariant under Eq. ( |3.1(J| ). We stress that from a formal point of view any (real) 
function represented by a power series would yield a realization of gauge invariance and 
that the choice of the electromagnetic form factor is entirely motivated on physical grounds 
through Eqs. (pD and |T|). 

Although a description of a finite-size pion based on Eq. ( |3.7| ) is mathematically con- 
sistent, it turns out to be inconvenient as soon as interactions of several particles with 
different form factors are considered. For instance, even a simple local interaction of fields 
(p a (a = 1, • • • , n) of the form 

O mt (<p a ) = gJl^, (3.14) 

a 

with a vanishing net charge associated with the product of fields, is no longer automatically 
gauge invariant as soon as different fields transform with different form factors F a . Under 
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the gauge transformation of Eq. (|3.10|) , the interaction Lagrangian, Eq. (|3.14|) , picks up a 



space-time-dependent phase factor 

Y[exp[-ieF a (-d 2 )A(x)} ± 1. (3.15) 

a 

On the other hand, this would not happen for the transformation law of Eq. ( |3.3|) , because 

l[exp[-ieZ a A(x)] = 1, (3.16) 

a 

provided the electric charge is conserved at the vertex, i.e., J2a Z a — 0. 
We therefore redefine the field as follows : 

<f>(x) = ip(x)exp[tef(-d 2 )d^(x)}, (3.17) 



whereQ 



f(q 2 ) = - 2 [F(q 2 )-Z] (3.18) 



generates a new field which transforms canonically under Eq. (|3.10| ), i.e. 

ip = exp[-ief(-d 2 )d^](j) (3.19) 
i — ^ exp[-tef(-d 2 )(d^ + d 2 A)] exp[-ieF(-d 2 )A](j) 
= exp[-ief(-d 2 )(d^ + d 2 A)} exp{-ie[Z - d 2 f{-d 2 )]A}<j) 
= exp[-ze/(-<9 2 )<9^] exp(-ieZA)(j> = exp(-ieZA)tp. (3.20) 

Rewritten in terms of (p, the Lagrangian (|3.13|) reads 

Co(Dfa <p) = Dfc{D f ^) ] - M 2 W \ (3.21) 

where the second (noncanonical) covariant derivative Dfa is defined as 

Dfo = [d, + ieZAp + ief(-d 2 )d»^ u ]<p, (3.22) 

and = d^A u — d u A^. This second definition ensures canonical gauge invariance and 
simultaneously accounts for the finite size of the particle. 

Organized in powers of the elementary charge, the Lagrangian can be expressed in terms 
of the canonical covariant derivative as 



6 Recall our assumption that F can be expanded in an absolutely convergent series. 



13 



£o(-D£y, <p) = ^o{D^ip, ip) + C^D^tp, <p) + C 2 {if). 



(3.23) 



C and 

Cxip^tp) = -ie^D^} f(-d 2 )d u ^ (3.24) 

generate the electromagnetic vertex 

F f i(p,p) = (p+p) f iF(q 2 )-q tl (p' 2 -p 2 )f(q 2 ), q = p - p, (3.25) 

which now satisfies the Ward-Fradkin-Takahashi identity 

q^(p',p) = Z{p' 2 -p 2 ) = Z[A-\p') - A-\p)] (3.26) 

with the free propagator. The last term 

= e 2 [f(-d 2 )d u ^][f(-d 2 )d^, P W (3-27) 

vanishes, when at least one of the photons is real, thus satisfying d u J rflu = 0. 

In the following we define the (generalized) Born terms of the virtual Compton scattering 
amplitude as the VCS amplitude constructed from either of the Lagrangians of Eqs. ( |3.7| ) 
and ( |3.21| ). According to the equivalence theorem of Lagrangian field theory [^,^], the 
scattering amplitude does not depend on the Lagrangian used as long as all external pions 
are on shell. To be specific, one obtains 



f$T = e 2 F(q 2 )F(q' 2 ) 



^ (2p + q),(2p' + q') u (2p - q%(2p' - q) t 



(p + q) 2 — M 2 (p- q') 2 - M 2 



(3.28) 



In fact, such a form of the Born amplitude is standard for a discussion of structure-dependent 
characteristics of the target, and we have shown that this is a very natural generalization of 
the Born amplitude for a point-like particle to the case of a finite-size particle. As discussed 
in Ref. |58] in detail, such a generalization incorporates all low-energy singularities of the 



total VCS amplitude so that the non-Born part of the amplitude can be expanded in powers 
of small photon momenta giving rise to (generalized) polarizabilities. 



IV. LOW-MOMENTUM EXPANSION AND GENERALIZED DIPOLE 

POLARIZABILITIES 



A well-known, general method of obtaining the low-energy expansion of a reaction am- 
plitude consists of expanding invariant amplitudes free from kinematical singularities and 
constraints in a power series (see e.g. p0| , p9| ). Following this method, we first decompose 
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the invariant amplitudes Bi[y 2 ,q ■ q',q 2 + q' 2 ,q 2 q' 2 ) into generalized Born and non-Born 
contributions, 

B i = Bf anx + Bf B , i = l,-..,5. (4.1) 

The generalized Born, or convection contribution is associated with a set of diagrams de- 
scribing single-pion exchanges in s- and w-channels with 77T7T vertices taken in the on-shell 
regime. As we have seen in the previous section, additional non-pole terms are necessary 
to render the generalized Born terms gauge invariant. The thus constructed amplitude pos- 
sesses all the symmetries of the total amplitude Tycs and contains all singularities of Tycs 
at low energies. Using the specific form of the Born amplitude Ty<2g n given in the previous 
section, we find the (generalized) Born parts of the invariant amplitudes Bi, 

B f°™ = (q ■ q')C, 5 2 Bom = -AC, C = . ^Ti^^L » ( 4 - 2 ) 



M 2 )(u-M 2 ) 



and Bf otn = for % = 3,4,5. 



At energies below inelastic thresholds, the non-Born parts of Bi are regular functions of 
the kinematical variables. They determine the deviation of Tycs from its Born value of Eq. 
In particular, when the momenta of both photons are small, q ~ q' — > 0, one obtains 

Tvcs = T^i + I W^O) + (P,Fn(P P F' pv )h(0) + 0(g 4 ), (4.3) 

where the constants 6<(0) = 5^(0,0,0,0), % = 1,2, can be related to the ordinary electric 
and magnetic dipole polarizabilities of low-energy real Compton scattering, 

8vrMa = -b x (0) - M 2 b 2 (0) , 8ttM/3 = bi (0) . (4.4) 



Equations (O) and ( |4.4|) provide a Lorentz-invariant form of the low-energy theorem for 



virtual Compton scattering up to second order in the photon momenta (for the case of the 
nucleon, see Ref. [pi). 

Now, following the original idea of Guichon et al. 0, we consider the case when the 
final photon is real and has a very small momentum q' 0, whereas the initial photon 
momentum q is allowed to be virtual and is not necessarily small. As may be seen from 
Eq. ( |2.29| ), the amplitude of 7*7r — > 77? with a real final photon is determined by three 
invariant amplitudes Bi, B2, and i?3,[] 



Tvcs = -zF^F'^Bi + (P,Fn(P p F' pu )B 2 + (PV^K^^J^. (4.5) 



"The definition of P of Eq. ( p.6[ ) differs by a factor of 1/2 from Ref. p7|]. 



15 



This equation has a particularly simple form in the pion Breit frame (PBF) defined by 
P = 0, i.e., p = —p', in which the initial and final pion are treated on a symmetrical 
footing. Introducing the Fourier components of the electric and magnetic fields as 

E = i(q e - qe Q ), B = iq x e, 
E' = -i(q' e'* - q'e'*), B' = —iq' x e'*, (4.6) 



we can rewrite Eq. (|4.5| ) in the Breit frame as 

Tvcs = [(B ■ B')B 1 - (E ■ E'^B, + P 2 B 2 ) + (E ■ q)(E' ■ q)P 2 B 3 ] pBp . (4.7) 

In general, even after subtraction of the singular (Born) parts of the amplitudes, all Bf* B 
still depend on the photon energies and the scattering angle and thus describe a series of 
multipoles and dispersion effects. However, in the limit of q' —>■ 0, only scalar structure 
functions depending on q 2 survive: 

b l (q 2 ) = Bf B (0,0,q 2 ,0). (4.8) 

These functions yield the non-Born parts of the coefficients in (|4.7|) which we interpret as 
generalized electric and magnetic dipole polarizabilities: 

8nM(3(q 2 ) = b 1 (q 2 ), 
8nMa T (q 2 ) = -h(q 2 ) - (V - ^ b 2 (q 2 ), 

^Ma L {q 2 ) = -hiq 2 ) - [m 2 - ^ [b 2 {q 2 ) + q 2 b 3 (q% (4.9) 

where P 2 has been taken in the limit q' = as well, i.e., P 2 = M 2 — g 2 /4. If the initial 
virtual photon has a transverse polarization in the Breit frame (E = E? -L q), the pion 
response to the transverse electric and magnetic fields is determined by the (generalized) 
transverse electric and magnetic polarizabilities: 

(T$? S ) T = 8ttM[(E t ■ E')a T (q 2 ) + (B ■ B'){3(q 2 )} + (higher orders in q'). (4.10) 

For a longitudinal polarization {E = El \\ q) in the Breit frame, 

(T^ C B S ) L = MM{E L ■ E')a L {q 2 ) + 0{q' 2 ). (4.11) 

In the real-photon limit, q 2 — > 0, we have 

0(0) = 0, a L (0) = a T (0) = a. (4.12) 

All thus defined polarizabilities are functions of q 2 free from kinematical singularities. In 
particular, for pion, kaon, or nucleon targets they are regular functions below the two-pion 
threshold, q 2 < Am 2 , and this region includes all space-like momenta. 
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We will now interpret a^q 2 ), ariq 2 ), and (3{q 2 ) by means of a semi-classical qualitative 
picture. To that end, let us consider a system of n polarizable (neutral) constituents at 
positions r a (a = 1, • • • , n) with electric polarizabilities a^. We allow the constituents to 
be anisotropic, i.e., the polarizability tensors are symmetric, afj = a% but not necessarily 
diagonal (oc 8ij). The system will respond to an external electric field E(t,r) by acquiring 
a dipole moment^ 

n 

Di(t) = 4tt (t,r a ). (4.13) 

0=1 

Slow oscillations of this dipole moment generate radiation of an outgoing long- wavelength 
electromagnetic wave E'(t, r) through the interaction — D ■ E. For an incoming plane wave 
with momentum q, i.e., E(r) = Eexp(iq ■ r), the amplitude for a transition to an outgoing 
plane wave with a very small momentum (viz. q'p <C 1, where p characterizes the system's 
extension) reads 

n 

f fi = 4vr £ a% exp(iq ■ p a ) , (4.14) 

0=1 

where p a = r a — R are the positions of the constituents with respect to the center of mass 
R of the system. The continuum limit of a system with spherical symmetry must be of the 
form 

n 

<j exp(iq ■ p a ) -> aij(q) = a^q)^ + a T (q){5 i:j - q^qj), (4.15) 

o=l 

where and do not depend on the direction q of q. In this way we recover the structure 
of the VCS amplitude given by Eqs. (pTTO| ) and flPip . 



If the system under consideration is exposed to a static and uniform external electric field 
E', an electric polarization "P is generated which is related to the density of the induced 
electric dipole moments: 

n 

Vi(r) =AtiJ2 <• 8 3 (r ~ r a ) E] = Ana i3 (r - R) £j . (4.16) 

o=l 

The tensor (Xij(r) is nothing else but the Fourier transform of the polarizability tensor of 

Eq. (wm n 



8 The factor of 4-7T is related with the (standard) use of Gaussian units for the polarizabilities but 
natural units for charges and fields. 

9 We do not use different symbols for a function f(t) and its Fourier transform f(uj). 
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aij{q) exp(-iq ■ r) 



dq 



(4.17) 



If we define 



B(q) 



«l(<?) - OL T (q) 



(4.18) 



Eq. ( 4.15 ) can be rewritten as 

atij(q) = a T (q)Sij + B(q)q i q j 



(4.19) 



such that the Fourier transformation results in 

atij(r) = a T (r)5ij - V»VjB(r). 



(4.20) 



Because of B — B(r), the second term of Eq. Q4.2U| ) reads 

- Vi V,B(r) = - (S„ - ^ - r -f±B"M. 



(4.21) 



On the other hand, from Eq. Q4.18 ) written in the form q 2 B{q) = ai{q) —atT(q), one obtains 
for the Fourier transform 



V 2 5(r) = -B"(r) - -B'Or) = a L {r) - a T {r). 

r 



(4.22) 



Eliminating B"(r) from Eq. ( |4.21 ) allows one to rewrite Eq. ( |4.20| ) as 



r,r,- 



a 



j(r) = a L (r)-^- + a T (r){8,, 



'J r 2 



+ 



(4.23) 



Finally, the last term of Eq. ( |4.23| ) is determined by reexpressing Eq. (|4.22 ) as 

d 



dr 



r B (r)] = — r [az,(r) — ar{r)]. 



(4.24) 



which, assuming the boundary condition \im. r _ ¥OQ r 2 B' (r) = 0, is solved byQ 



10 Instead of Eq. Q4.25p we could also use 



r 2 B'{r) 



r [tti(r') — «T(r')] dr' , 



resulting from the boundary condition lim^o r B'(r) = 0. Both results are identical, because 

4irr 2 [aL(r) — aj-(r)] dr = ai(q = 0) — ax{q = 0) = 0, 



where we made use of Eq. (gig ). 



r 2 B'(r) = r' 2 [a L (r')-a T (r')}dr'. (4.25) 

Jr 

In other words, given the Fourier transforms 

«L,r(r) = J a L , T (q) exp(-iq ■ r) (4.26) 

the density of the full electric polarizability of the system, = J2 a a ij> can be reconstructed 
as 



3f*-f* • ■ ■ poo 

aij (r) = a L (r)rirj + a T (r)(5 tJ - r^) + - 3 l] j [a L (r') - a T (r')] r' 2 dr'. (4.27) 



In this context, it is important to realize that both longitudinal and transverse polar- 
izabilities, ai and oct, respectively, are needed to fully recover the electric polarization V '. 
The longitudinal polarizability is special, though, because it completely specifies the induced 
polarization charge density of the system, 

Sp( r ) = -V • V(r) = -Att(E' ■ V)a L (r), (4.28) 

where we made use of 

Viay(r) = Vja L (r) (4.29) 

which follows from Eqs. (|4.15|) and ( [4.17| ). Combining partial integration with the divergence 
theorem, one finds for the Fourier transform of the induced polarization charge 

6p(q) = J 6p(r) exp(iq ■ r) dr = iq ■ V(q) = Ama L (q) q • E' . (4.30) 

Such an induced charge density is the source of a longitudinal electric (Coulomb) field and 
thus generates an effective coupling of the type El ■ E' . 

At the same time, the transverse polarizability olt describes rotational displacements of 
charges which do not contribute to Sp(r). They can generate electric fields only for a finite 
frequency, q' Q ^ 0. Therefore, in the limit q' —> the effective coupling Ex ■ E' should vanish 
faster than El ■ E'. This is indeed the case, as we will see below. 

The relation ( |4.30|) suggests an intimate connection between the longitudinal polarizabil- 
ity ai and the charge- density operator of the system. In a forthcoming publication we will 
verify this in a quantum-mechanical and fully relativistic framework. 

Similar considerations apply to the magnetic part of the VCS amplitude. In this case the 
electric polarizabilities should be replaced by the magnetic ones, Since the magnetic 
induction is always transverse (i.e., B q = 0), terms containing q^j in the magnetic analogue 
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of Eq. (|4.15|) do not enter any observable and can thus be omitted. Hence, the unobservable 
"longitudinal" magnetic polarizability (3l{q) can for a h Q be chosen to be identical with the 
transverse one, flr{q) rather than only at the point q = 0, where the equality /?l(0) = /Jr(0) 
is dictated by analyticity. With this choice /%(<?) = /3(q) and the analogue of Eq. (|4.27|) 
reads 

fair) = PWn. (4.31) 

Then the magnetization induced by the uniform external magnetic field and the corre- 
sponding induced electric current Sj(r) = V x At(r*) are, up to an arbitrary gradient, given 
in terms of the density of the magnetic polarizability as 

M{r) = 47r(3(r - R) B' , (4.32) 

where the Fourier transform of /3(r) is nothing but the generalized magnetic polarizability 

/%): 

/3(r) = J /3(q) exp{-iq • r) ^L. (4.33) 



Let us conclude this section by recalling that Eq. ( |4. 10| ) was obtained by keeping the 



lowest multipolarities of the final (soft) photon. This, however, leads to different powers of 
the photon momentum q' in such an expansion. For q' — > 0, the energy of the initial photon 
in the Breit frame vanishes as well, because go — q'o in that frame. It then follows from 
Eq. ( |4.6| ) that the transverse electric field Et is of higher order in q' than the (transverse) 
magnetic field B. This is also clear from the relation q 2 Ex = —qoQ x B. When only terms 
up to order 0(q') are retained, the transverse electric field does not contribute! In order to 
translate Eq. ( |4.7|) into a power expansion of the non-Born part of the VCS amplitude up 
to 0(q' 2 ), one has to add two more terms proportional to (q ■ q') [dBf B / d(q ■ q')] q >=o an d to a 
similar derivative of the function Bi + P 2 B 2 . These terms introduce an additional angular 
dependence and therefore higher multipoles (quadrupoles). 

The physical process of photon electroproduction e(k) + h — > e'(k') + h' + j(q') consists 
of the Bethe-Heitler contribution, in which the real final photon is emitted by the initial 
and final electrons, respectively, and the VCS contribution. The virtual photon of the VCS 
part, 7(g), is described in terms of the polarization vector 

g 

e M = -^u e (k')^^u e (k), q = k-k', (4.34) 



which is determined by the electron-scattering kinematics. For q' — * 0, such an remains 
finite. Therefore, the transverse electric field Et = iq^T created by the electron transition 
current is of order 0(q') in the Breit frame and is suppressed in comparison with the magnetic 
and longitudinal electric fields generated by the current. As an immediate consequence 
the non-Born part of the VCS amplitude to order 0(q') is characterized by two structure 



20 



functions [viz. ai{q 2 ) and (3{q 2 )} rather than by all three functions appearing in the dipole 
approximation. 



Although we arrived at this conclusion by considering the VCS amplitude in the Breit 
frame [p = (q' — q)/2], it is clear that two independent structure functions characterize the 
amplitude to order 0(q') in any other frame such as, for example, the center-of-mass (cm.) 
frame. This is true because the amplitude itself is Lorentz invariant and at the same time 
the real-photon energy, uo 1 = \q'\, remains of the same order 0(q') for any finite Lorentz 
boost. 

The above consideration gives a transparent explanation of the theorem established in 
Ref . , which states that there are only two independent structure functions which deter- 
mine 0(q') terms in the so-called full amplitude T^ cs of virtual Compton scattering for a 
spin-0 particle in the cm. frame. Here, "full" refers to the fact that the polarization and 
the intensity of the initial photon are given by Eq. (|4.34|) . 

In terms of the notation introduced by Guichon et al. ||, this theorem establishes a 
linear relation between the cm. generalized polarizabilities p(° 1 ' 01 )° ) p(n>ii)o^ anc j p(oi,i)o^ 
leaving only two of them independent. In the cm. frame, the transverse electric field is not 
vanishing, and the above theorem can be re-phrased as establishing a linear combination of 
the El, Et and B responses in the cm. frame which vanishes when q' = 0. This is just the 
Et response in the Breit frame. 

The explicit relations between the cm. polarizabilities and the quantities ai{q 2 ) and 
(3{q 2 ) appearing at order 0(q') read 



where E cm = M-q 2 /{2M) and q cm = ^-q 2 + q A j{^M 2 ) denote the energy and the absolute 
value of the three-momentum of the initial pion in the cm. frame at threshold (</ = 0). The 
spin factor 2J + 1 removes a related spin dependence hidden in the quantities P^ p L ,pL ^ s and 
is needed when our "pion" represents a spin-averaged hadron of spin J ^ 0. 



When considering the Fourier transforms, the additional factor of \fE cm in Eq. ft4.35|) and 



the use of the Breit-frame momentum transfer gereit — V — <? 2 instead of the cm. momentum 
transfer q cm will generate a difference for the spatial distributions, especially for such a light 
particle as the pion. From the analogy with the well-known case of electromagnetic form 
factors, where spatial distributions are obtained using the Breit-frame variables, we expect 
that a meaningful Fourier transformation in the case of the generalized polarizabilities also 
requires the Breit frame. Indeed, the only difference between the kinematics of the reaction 
of VCS, 7*7r — > 7 ; 7r, and the kinematics of the reaction 7*7r — > tt, in which the form factors 
are studied, originates in the presence of an additional photon 7' which carries a vanishing 
momentum q' = 0. 
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In analogy to Eq. ( |2.29| ), the structure-dependent effects as seen in VCS with one (soft 
or hard) space-like virtual and one soft real photon can be encoded in the following effective 
Lagrangian: 



£poia riz = fxT^T^H + 2 [kF a ^ - 2b 3 ^ dpffTf^PaPtyfa (4.36) 

where hi are differential operators acting on the electromagnetic fields which are determined 
by their Fourier components h(q 2 ) [see Eq. ( |4.9|) ]. The above Lagrangian contains all possible 
gauge-invariant terms involving at least one field strength tensor without derivatives 
[0(u% 



V. SPATIAL DISTRIBUTIONS IN CHIRAL PERTURBATION THEORY 

A. Preliminary remarks 

In accordance with the ideas presented in the previous section, we now consider the 
Fourier transforms of the g-dependent polarizabilities, generically denoted by F(q 2 ), and 
discuss the corresponding spatial distributions F(r). 

There is a well-known objection against a straightforward interpretation of such F(r) as 
a spatial distribution. The argument is related to the fact that the velocities of the target 
before and after the interaction with the virtual photon depend on the photon momentum. 
If we think of the target as a composite system of, say, quarks we would expect that the 
matrix element 

J ip}(R' \t')(R! 'y\0\R,T)iJ;i(R,T) dRdr dR' dr (5.1) 

of a transition operator like O — j j fl (x)A fl (x) d 3 x for a one-photon reaction^ depends on 
both the internal (Lorentz-invariant) variables r of the pion and on the pion's center-of-mass 
variable R (cf. Ref. [(hj). Since a relativistic wave function ip(R,r), in general, does not 
factorize into a product of the type 0(r)exp(ip • R), where 0(r) denotes a p-independent 
internal wave function, some part of the full g-dependence of the transition matrix element 
may be related with the cm., partly kinematical effects of p on 4>{r). We come closest to 
associating the Fourier transform F(r) of the full matrix element ( |5.1| ) with the internal 
spatial structure of the particle by evaluating this matrix element in the Breit frame, in 
which the pion is at rest on the average, i.e., p + p' = 0. 



11 Another example would be O = JJ j u (x')A' u (x')G e(x' , x)j >l (x)A fl (x) d 3 x d 3 x' involving two cur- 
rents, where Ge is the Green function. 
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There is a simple phenomenological argument suggesting that the generalized magnetic 
polarizability (3{q 2 ) defined according to Eq. Q4.9Q is indeed only related to the internal 
structure. The point is that this quantity is a function of q 2 without kinematical singularities. 
In other words, any irregularity in its q 2 behavior is not caused by a Lorentz contraction 
and has nothing to do with singular 7 factors or with quantities like Po — \\f^M' 2 + q 2 . 
Moreover, even singularity- free quantities like Pq = P 2 should be irrelevant, because the 
momentum scale on which the amplitude b\(q 2 ) changes has nothing to do with the particle 
mass M itself and is fully determined by interactions. 



For the other polarizabilities the situation may be more complicated. The sums a^(g 2 ) + 
P(q 2 ) and ar{q 2 ) + P(q 2 ) contain an overall factor of P 2 [see Eqs. Qi.9|) l which, for the 



pion, introduces a small kinematical mass scale into the q 2 behavior of these sums. So, 
perhaps a more meaningful consideration of spatial distributions relating to a^(r) + /3(r) 
and ax(r) + (3{r) is obtained with the factor P 2 /M 2 excluded from a^q 2 ) + (3{q 2 ) and 
otriq 2 ) + fi(q 2 ) before performing the Fourier transformations. In the following discussion 
we will not encounter this problem, because in the theory considered, namely, ChPT at 
lowest nontrivial order, either the sums otL^q 2 ) + (3(q 2 ) and ar{q 2 ) + P(q 2 ) are exactly zero 
[for pions and kaons at 0(p 4 )} or the particle mass M is considered to be infinite [for baryons 
in HBChPT at 0(p 3 )}. 

We will illustrate, by means of the more familiar example of form factors, that associating 
a generic F(r) with the internal structure of the particle leads to a self-consistent picture 
and does not create any visible problems even in the case of such a light particle as the 
pion, for which the relativistic interlace of cm. and internal variables is maximal. To be 
specific, we will consider form factors calculated in the framework of ChPT, mainly for two 
reasons. First, we want to discuss the generalized polarizabilities F(q 2 ) predicted by the 
same theory in order to check that our consideration of polarizabilities at scales r ~ l/m n 
is not in conflict with other observables. Second, at present ChPT is the best tool for 
describing hadron structure at large scales and it is the only theory which agrees with the 
recent MAMI data on generalized polarizabilities of the proton [TT| . 



We would like to mention the following technical aspect concerning the Fourier trans- 
formation of a q distribution obtained within ChPT. Such distributions are only reliably 
known for small momenta, q = 0(m n ). Moreover, a straightforward integration over q does 
not exist, because the integrand typically diverges for large values of q. We therefore enforce 
convergence by introducing a cutoff A. Clearly, such a cutoff disturbs the corresponding r 
distributions at distances r < 1/A which are beyond the scope of ChPT. On the other hand, 
one might expect that the results are cutoff independent when r ^> 1/A. In practice, we 
make use of a Gaussian cutoff and, for any g-dependent form factor or polarizability F(q 2 ), 
we calculate F(r) as 

rv x y A [°° P / n2 x sin(Qr) , Q 2 ^Q 2 dQ 

F(r)=hm 4^ F(-Q ) exp (--) _. (5.2) 

Depending on how small r is, we have to choose A large enough in order to approach the limit 
of A = 00. In particular, for all generalized polarizabilities considered below we have found 
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the regularized Fourier integral of Eq. ( |5.2|) to be independent of A for A > 30 GeV within 
an accuracy better than 2% even at distances as short as r = 0.1 fm. Stated differently, 
the cutoff A = 30 GeV is sufficient to resolve spatial distributions of polarizations at scales 
r ~ 0.1 fm. In the case of electromagnetic form factors, having steeper spatial distributions 
(see below), the 30 GeV cutoff is sufficient for a good resolution up to distances of r ~ 0.2 
fm. 



There is yet another way of calculating the Fourier integral of Eq. (|5.2j ) based on a 
contour deformation in the complex Q plane. This method is applicable when the analytical 
continuation of F(q 2 ) to time-like momenta is known as in the case of the ChPT predictions. 
Since all singularities of F(t) are located at real positive t, it is possible to write a dispersion 
relation for FjQ 

F(?2) = ?£ ImF(t) I^bF < 5 - 3) 



which allows one to recast the Fourier integral for F(r) at r > as a superposition of Yukawa 
functions: 

1 f°° r 

F(r) = — \ e~ rVl ImF(t) dt. (5.4) 

'-±71 T J L; n 



We made use of both methods and arrived at identical results for F(r). 

It is worthwhile recalling that polynomial pieces in F(q 2 ) generate 5(r) terms or deriva- 
tives thereof in the Fourier transform F(r). Such terms typically originate from higher-order 
counter terms in the Lagrangian which are multiplied by a priori unknown low-energy con- 
stants. However, in the Fourier transform, they do not contribute to F(r) at finite r ^ 0. 
In other words, the Fourier integral acts as a filter which only transmits genuine effects of 
pion (or kaon) loops through their contributions to a nonpolynomial part of F(q 2 ) and to 
ImF(g 2 ), respectively. 



B. Form factors 



As a first illustration, we briefly discuss the scalar and vector form factors of the pion 
as obtained in two-flavor ChPT in the limit of isospin symmetry. These form factors 
parametrize matrix elements of the scalar density S(x) = rh[u(x)u(x) + d(x)d(x)] with 
m = m u = m d and the isovector electromagnetic current jY(x) = |<?( x ) r 37 nQ{ x ) > respec- 
tively: 

(n(p')\S(0)Hp)}=F s ( q 2 ), 
(K + (P'M(0)W + (P)) = (p + p')„Mq 2 ), q = p'~p. (5.5) 



Additional subtractions may be required resulting in additional polynomial contributions. 
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Recently, the one- loop calculations of Fs(q 2 ) and Fy(q 2 ) by Gasser and Leutwyler have 



been extended to the two-loop level 0{p ) [|42| , f43 



In order to simplify the discussion, we perform two subtractions in the form factors F(t) 
and plot subtracted (and normalized) functions, 



F(0) 



F{t)-tF'{0) - ^t 2 F"{0) 



(5.6) 



By that means we avoid polynomial contributions of 0{p i ) and 0(p 6 ), which depend on 
low-energy constants, and emphasize the pieces originating from pion loops.0 As stated 
before, such subtractions are not visible through the Fourier filter at r > and are thus 
irrelevant for the determination of F(r). 



At 0(p 4 ), the subtracted scalar form factor of the pion reads 



Fs(q 2 



ml 



lQn 2 F2 



2x 



1 _ f(Vl/ . 19x 2 
-J {0 >{x) + - 



lOx 



120 



(5.7) 



where x = q 2 /m 2 , and the function J(°>(x) is defined as^ 



j(°) ( x ) = / lnfl + x(y 2 - y) - i0 + ] dy = - 2 - a In 
Jo 



a- 1 



a = Wl 



4 

j 

x 



x < 0. 
(5.8) 



As numerical values, we use F n = 92.4 MeV [4~7] and the charged pion mass = 139.6 
MeV. The polynomial in Eq. ( |5.7| ) results in vanishing first and second derivatives of F$ at 
t = 0. 



At one-loop order, the subtracted vector form factor is of a similar form: 

^ fx -4 , 3x 2 -20x\ 



Fv(q 2 



1 



ml 



167r 2 F 2 V 6 



180 



(5.9) 



At the two-loop level, the scalar and vector form factors are given by more lengthy expres- 



sions which can be found in Refs. |42|,43|. To be specific, we made use of Eqs. (3.6)-(3.8) and 
(3.16)-(3.18) of Ref. [Q, using the parameters (low-energy constants) Ti = —1.7, l 2 = 6.1, 
l 3 = 2.9, h = 4.472, r 6 = 16.0 (set I @). 



13 Alternatively, we could keep the polynomial contribution of the pion loops. However, in that 
case the result would still depend on the renormalization condition chosen. 



14 



The results for < x < 4 and 4 < x are obtained by analytical continuation. 
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For comparison, we discuss as another example the isovector electric form factor of the 
nucleon, G^(q 2 ), to leading nontrivial order [C(p 3 )] within HBChPT With the above 

two subtractions one obtains^ 

<9\ ^-8 j(0)(y) + 21^-40x \ (51Q) 



16tt 2 F2 V 6 v ' 180 

where = 1.267 is the axial-vector coupling constant. 

The scalar form factor of the nucleon and the corresponding spatial distribution were 
recently discussed by Robilotta [§|| in the context of the two-pion-exchange contribution to 
the NN potential (see Fig. 8 of that reference). 



In Fig. [I] we show the three form factors F(q 2 ) of Eqs. ( |5.7|) , (|5.9|) , and (|5.10|) and their 



corresponding Fourier transforms F(r) which, in the case of F — Fy and F = G E , are 
interpreted as the electric charge density.^ For the sake of comparison we also present the 
results of a pole approximation to those form factors, 

2 

777 

F(q 2 ) = — Anr 2 F(r) = m 2 r e~ mr ', (5.11) 

m 2 — q 2 

with the mass scale m = 770 MeV arbitrarily chosen to be the p-meson mass. 

As a first observation, we note that all spatial distributions considered are positive at 
the one-loop order. Such a positive sign confirms a (maybe too) naive understanding of the 
peripheral structure of the target as being created by a cloud of virtual pions which create, 
in the case of a tt + or proton target, a positive electric charge density of finite size. Also the 
scalar density seen by an external scalar field carries a positive sign. The QCD coupling of 
an external scalar field s(x) to the scalar density S(x), 

£ ex t = —rh[u(x)u(x) + d(x)d(x)]s(x), (5-12) 

is described through the Lagrangian ^-F 2 Tr (xU^ + Ux^) at lowest order in ChPT, where, 
in the present case, x = m 2 s(x). Inserting U = (er + it ■ ir)/F, the lowest-order interaction 
reads 

££ = -~mln 2 (x)s(x) = -S cS (x)s(x). (5.13) 



15 In this case one subtraction would be sufficient to remove low-energy constants. Note that the 
functions J(q 2 ) and C(l 2 ) used in Ref. [Q are related with the function j(°\x), Eq. ( |5.8|) , by 
167r 2 C(g 2 ) = -J {0) (x) and (96vr 2 /m 2 )J(g 2 ) = -2x/3 + (x - 4)j(°)(x). 

16 Analytical representations of F(r) in ChPT at one-loop order are given in Appendix [A]. 
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0.2 0.4 0.6 1 2 

-q (GeV 2 ) r (fm) 



FIG. 1. Left panels: the (twice-subtracted) scalar and vector form factors of the pion in ChPT 
at order 0(p & ) |^3| and the isovector charge form factor of the nucleon in HBChPT at order 
C(p 3 ) p4| , f45| ]. Right panels: corresponding spatial distributions obtained as Fourier transforms. 
Dashed lines: one-loop prediction. Solid lines: two-loop prediction (for pions only). Dotted lines: 
prediction of the pole dominance model due to a scalar or vector meson of mass 770 MeV. 
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Thus, S eS (x) has manifestly positive matrix elements which, when probed through s as part 
of the pion cloud, lead to a positive density. 

When two- loop corrections are taken into account, the spatial distributions at small 
distances, r < 1 fm, change drastically. In particular, the charge density Fy{r) of the pion 
which at the one- loop level was concentrated at r < 0.3 fm, now extends up to r ~ 0.6 
fm. An even more dramatic effect is observed for the scalar density Fs(r), which due to the 
two-loop contribution changes sign at r < 0.7 fm. However, one should keep in mind that 
ChPT is a low-energy effective field theory and it is likely that higher-order corrections will 
change the spatial distributions at short distances substantially. An indication for such a 
scenario is given by the dotted lines in Fig. |l] which refer to short-range mechanisms such 
as vector mesons. Of course, such mechanisms are only included in ChPT by means of 
low-energy couplings in the effective Lagrangian. 

Finally, let us emphasize again that the above examples do not pretend to prove that 
spatial distributions F(r) have a strict operational sense. However, they certainly suggest a 
simple intuitive picture of the hadron periphery which is the natural domain of ChPT. 



C. Polarizabilities 



We now extend the discussion to spatial distributions associated with the generalized 
polarizabilities as introduced in Sec. [TV]. The predictions for the generalized dipole polar- 
izabilities cti{q 2 ) and j3{q 2 ) of the nucleon in the framework of SU(2)j HBChPT at 0(p 3 ) 
read IKIL 



a L (q 2 ) = E 2 W ['[Ss 2 + sml(lQy 2 - lQy + 9) - Zm\{2y - 1) 
Jo 



dy 



32s 5 /2 



E 



7T 



^■m n a(a + 1) 



arctan ^fa 



[a + l)(2a - 1) + 2a + 1 



a 



(3(q 2 ) = El f\2y - l) 2 (s + ml)(As - 3m 2 n ) 
Jo 



dy 



l Qs 5/2 



E 2 r, w n arctan Ja i 

(a + l)(2a + l) 2a- 1 , (5.14) 



16m T a(a + 1) L ^/a 

where s = m 2 r + q 2 (y 2 — y), a = —q 2 /(4:m 2 ), and E w = e^\/2/(87ri^) is the Kroll-Ruderman 
amplitude of 7r =l= photoproduction in the chiral limit .0 



17 The integral representations of Eq. ( 5.14|) are, of course, less convenient than the equivalent 



elementary formulas. They are given here only as a historical reference, in the form in which 
they were first calculated and reported |4q| . The same results were found independently by the 
authors of Refs. [p3|-21|. Recently, results have also been given in the framework of the Small Scale 



Expansion including the A-isobar as a dynamic degree of freedom [^] . 
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The transverse electric polarizability of the nucleon remains yet to be determined. 

There are no published calculations of the nucleon's generalized polarizabilities ai{q 2 ) 
and f3(q 2 ) within SU(3)j HBChPT, except for the case q 2 = considered in Refs. |49| , |50| . 
However, it is straightforward to extend the results of Eq. ( |5.14j ) to the SU(3)/ case in the 
limit of equal N, A, and E masses. In this limit, the square baryon-octet mass differences 
is considered as small in comparison with the square kaon mass, m K , which empirically is 
a good approximation. Then, the structure of Feynman diagrams contributing to Compton 
scattering in the SU(3) / and SU(2) / cases, respectively, is very similar, the only difference re- 
sulting (a) from different meson-baryon couplings for charged mesons (or, stated differently, 
from different Kroll-Ruderman amplitudes) and (b) from different masses of the Goldstone 
bosons entering the appropriate loop diagrams. The generalized polarizabilities of the nu- 
cleon in SU(3)/ HBChPT receive, in addition to the results of Eqs. Q5.14Q , contributions due 
to kaon loops which are given by the same expressions as in Eqs. ( p. 14j ) after the replacements 
m n —>■ rriK and E n —> Er, where 

Et = E 2 KR ( lN ^n±N) = (^) 2 L, 

E 2 K = E^N - K+A) + El^N K+Y.) = {^fltc (5.15) 



and 



t 2 m i p\2 t \\D 2 + 2F 2 , proton . , 

U = g A = (D + F), Ik = U d _ f)2 ^ neutron (5.16) 



(cf. ||49| , |50|| ) . In an 0(p 3 ) calculation, the difference between the pion and kaon decay con- 
stants is of higher order. Empirically, Fr — 1.22F n [47], but in our numerical analysis we 
make use of a universal value Fx = F n with F n = 92.4 MeV. Furthermore, we insert the 
SU(3)/ ratio F/D = 2/3. 

By applying the same procedure to the other octet baryons, one obtains a generic polar- 
izability F(q 2 ) as 



F(q 2 ) = Eq. (5.14) 



x 



9\ 



+ 



Eq. ( gng ) 



Ik 

x — . 



(5.17) 



The SU(3)/ coefficients I n and Ik are identical with the ones given in Eq. (4) of Ref. j49| for 
the case q 2 = 0. For convenience, we collect these coefficients in Table | (using F/D = 2/3). 
Note, however, that the pion-loop contribution for A and £ is perhaps not given very 
accurately by this procedure, since the baryon mass difference in the transitions 7r ± A <-> £ 
is not fully negligible in comparison with the pion mass m n . 

At 0{p A ), SU(2)j and SU(3)/ ChPT predictions for all three dipole polarizabilities 
a L{q 2 ), c^riq 2 ), and (3(q 2 ) have been reported for pions and kaons [p7p8| . At that order, 
the three generalized dipole polarizabilities are degenerate, i.e. 
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TABLE I. Flavor coefficients 1^ and Ik determining pion and kaon loop contributions to the 
generalized polarizabilities of octet baryons in SU(3)j HBChPT at order 0(p 3 ). 



B I v 


In/g 2 A 


Ik 


Ik/9 2 a 


V (D + F) 2 


1.00 


Id 2 + 2F 2 


0.56 


n {D + F) 2 


1.00 


(D-Ff 


0.04 


A §L> 2 


0.48 


\D 2 + 3F 2 


0.60 


£+ ID 2 + 2F 2 


0.56 


\d + f) 2 


1.00 


S° 4F 2 


0.64 


D 2 + F 2 


0.52 


E- ID 2 + 2F 2 


0.56 


(D - F) 2 


0.04 


E° (D - F) 2 


0.04 


(D + F) 2 


1.00 


E' (D - F) 2 


0.04 


ID 2 + 2F 2 


0.56 


a L (q 2 ) = a T (q 2 ) = 


-fl(q 2 ) and ai(r) = ax( 


r) = -P(r), 


(5.18) 


and hence we only need to discuss, say, the generalized ma 


gnetic polarizability which can 


be expressed as 28] 








2 e 2 1 
™ ' ~ AttM (4ttF t ) 2 


A + (B w + J {0) '(x n ) + C K x K J {0) '(x K ) 


(5.19) 


Here, M is the mass of the particle in question, 








q 2 q 2 
m K 




(5.20) 


the function (x) is given by 








jV)'(x) 


_ J<®(x) _ 2J(°)(x)+x 
dx x(x — 4) 




(5.21) 


and A, B n , CV, and Ck are constants given below. The terms in Eq. (pM9j) depending on 
and xk represent contributions of pion and kaon loops, respectively. The g 2 -independent 


term proportional to A originates from a contribution at short distances represented by 


low-energy constants in the effective chiral Lagrangian 51 : 







A(tt ± ) = A(K ± ) = 64n 2 (L r g + L r 10 ) = — - = 0.90 ± 0.12, 

Fy 

A(tt°) = A(K°) = A(K°) = 0. (5.22) 



The numerical value of A(n ± ) is fixed |52j by the experimentally known axial (Fa) and 
theoretically known vector (Fy) form factors of the radiative pion decay 7r + — > e + v e ^j using 
F A /F V = 0.448 ± 0.062 @. The relation A(n) = A(K) is valid in the SU(3)/-symmetry 
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FIG. 2. Generalized magnetic polarizability (3{q 2 ) of pions and kaons at C(p 4 ) [^J2|. Dashed 
lines: contribution of pion loops. Solid lines: contribution of pion and kaon loops. Dotted lines: 
VMD predictions normalized to (3 as given by SU(3)j ChPT. 



limit. The other constants entering Eq. ( |5.19| ) are flavor-dependent coefficients which deter- 
mine contributions of pion and kaon loops, respectively: 



-1 for 7T°, 

(X I -1/2 for 7r± 

-1/4 for kaons, 



for 7T°, 

otherwise, 



C 



K 



-1/2 for AT*, 
—1/4 otherwise. 



(5.23) 



The generalized polarizability of the Goldstone bosons are shown together with the results 
for the proton, the S~, and the S~ in Figs. and |3|. 

The spatial distributions calculated as the Fourier transforms of the generalized polariz- 
abilities are shown in Figs. f| and |5|. The corresponding analytical results obtained through 
Eq. Q5.4Q are given in Appendix For large r, the pion and kaon loop contributions to a 
generic polarizability F(r) follow an exponential behavior e~ 2m ^ r and e~ 2mKT , respectively, 
as determined by the nearest singularities of F(q 2 ) for time-like momenta, q 2 = 4m 2 and 
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FIG. 3. Generalized longitudinal electric and magnetic polarizabilities a^(g 2 ) and (3(q 2 ) of the 
proton, the £~, and the S~ at 0(p 3 ) (see [21,48| and the text). Dashed lines: contribution of pion 



loops. Solid lines: contribution of pion and kaon loops. Dotted lines: VMD predictions normalized 
to a and j3 as given by SU(3)/ ChPT. 
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FIG. 4. Density of the magnetic polarizability /3{r) of pions and kaons at 0(p 4 ). The 



(^-singularity at r = is not shown (see the discussion of Eq. (5.24) in the text). Dashed lines: 
contribution of pion loops. Solid lines: contribution of pion and kaon loops. Dotted lines: VMD 
predictions normalized to (3 as given by SU(3) f ChPT. 
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FIG. 5. Density of the longitudinal electric and magnetic polarizabilities ai{r) and /3(r) of 
the proton, the S~, and the E~ at 0(p 3 ). Dashed lines: contribution of pion loops. Solid lines: 
contribution of pion and kaon loops. Dotted lines: VMD predictions normalized to a and (3 as 
given by SU(3)/ ChPT. 
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q 2 = Am 2 K . The 5-singularity at r = cannot be seen in these plots. However, at least for 
mesons, such a singularity exists for sure within ChPT, and it is determined by the asymp- 
totic value of the polarizability for q 2 — > — oo. Thus, the integrals of the spatial distributions 
over r > are 

lim / Anr 2 (3(r) dr = f3(q 2 = 0) - /3{q 2 = -oo). (5.24) 

£-►0+ J e 



The generalized polarizabilities of the octet baryons, given by Eqs. ( |5.14j ) and ( |5.19 ), vanish 



at infinity, so that the integral ( |5.24 ) gives just (3 in this case. This is not the case for 



mesons, since xJ^'{x) — > 1 for x — > — oo. The corresponding 5-singularity in 47rr 2 /3(r) is 
driven by the limits |53| 



lim /3 7r ±(g 2 ) 

q l ^> — oo 


9 - 


lim (3 w o(q 2 ) 


15 R 


lim (3 K ±(q 2 ) 

q^^ — oo 


R l 9 mn R 
2 rriK 


lim (3 K o(q 2 ) 

q z ^—oa 


= 3 — —(3 n o. 
m K 



(5.25) 



The SU(2) / results for pions are found in Eqs. (39) and (40) of Ref. [^7j . 



It is very interesting that the loop contributions behave exactly as one would expect from 
a classical interpretation of the Langevin diamagnetism. In such a picture, a change in the 
external magnetic field would produce circulating currents induced in the charge density of 
the meson cloud which on their part give rise to an induced magnetization. Both pion and 
kaon clouds are seen to be diamagnetic, at least at large distances r > l/m w . Simultaneously, 
these clouds generate a positive sign for the electric polarizability. However, in the case of 
the nucleon the diamagnetic character of the pion cloud disappears at distances r < l/m n , 
where paramagnetism prevails and j3(r) becomes positive. 

When a hadron is probed by photons of very small (space-like) momenta \q\ m n , the 
magnetic response is essentially only sensitive to (3{q 2 = 0) = (3. When the momentum \q\ 
increases, the very peripheral, negative part of the magnetic polarizability of the pion cloud 
no longer contributes to 

p(q 2 ) = J (3( r ) exp(iq ■ r) dr, (5.26) 

due to the oscillatory behavior of the integrand for large distances. This explains why the 
magnetic polarizability (3{q 2 ) for all the particles (pions, kaons, nucleons, etc.) universally 
gets a positive increase (see Figs. ^ and |3|). The slope of (3{q 2 ) at the origin is proportional 
to the mean square radius of the spatial distribution, 

d/3(0) 2tt r» N , . 
-^ = -1 ^(r)*. (5.27) 
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Obviously, for all mesons, the cloud distribution 4nr 2 3(r) < for all r > 0. Hence, the slope 
is negative as a function of q 2 (positive, when plotted against — q 2 ). Also, the curvature as a 
function of q 2 is concave. For the proton, the £~, and the S~ the asymptotic negative pion 



and kaon tails in the integral ( 5.27 ) dominate over the positive contribution coming from 
distances smaller than 1 fm, 1 fm, and 0.4 fm, respectively. This makes the slopes of the 
magnetic polarizability of all baryons considered positive as a function of —q 2 as well. On the 
other hand, such behavior is in some cases opposite to that expected from VMD. Suppose 
the q 2 dependence of 3(q 2 ) was determined by the porw mesons mediating electromagnetic 
interactions, 

2 

f3(q 2 )} =-^^— (5.28) 

r ^ J VMD y m 2 _ q 2 V ) 

From this, the spatial distribution of 3{r) 

A7rr 2 (3(r)} = m 2 r e~ mpr 3 (5.29) 

v J VMD P r v ' 

would have the same the sign as f3 for all r > 0. For the 7r° and for the octet baryons, 
having positive j3, such a sign is in conflict with the ChPT behavior [see Figs. |3| and |5| in 
which the VMD distributions are also shown]. From a phenomenological point of view, the 
full magnetic polarizability 3 of the neutral pion should be approximately three times the 
0(p 4 ) prediction, mainly due to a paramagnetic contribution of the Ml transitions 7r° — > uj 
or fp. In other words, a more realistic VMD curve should be three times higher than that 
shown in Fig. |], so that the above conflict with ChPT would be even more severe. 

In the case of the (longitudinal) electric polarizability of octet baryons the long-range, 
peripheral part is suppressed with increasing — q 2 . This part is relatively large for all baryons 
except for the S's, and thus aiiq 2 ) shows a rapid decrease with increasing space-like q which 
is steeper than for VMD [see Fig. || . 

Another instructive feature of the plots shown in Figs. (| and |^ is that the kaon-loop 
contribution is concentrated at about the same scale of r (or even at smaller scales) as 
the VMD contribution due to the large mass of kaon pairs, 2Mk ~ 1 GeV. Kaon-loop 
contributions exactly fall into the short-range region where one finds other contributions of 
similar range treated via low-energy constants of the effective Lagrangian of ChPT. However, 
at 0(p 3 ) no such counter terms contribute to the generalized polarizabilities. In other words, 
even though the relevant counter terms are formally of higher order in the power counting of 
ChPT, one may wonder whether their quantitative importance is underrated as compared 
to the kaon-loop contributions. Therefore, quantitative conclusions drawn from calculations 
keeping kaon loops and ignoring short-range contributions have to be treated with some 
care. For a similar conclusion, see Ref. f54 |. 



The spatial distribution of 3(r) for the nucleon shown in Fig. [5] may also shed light on 
an old question which has remained open for more than 30 years: Why is the magnetic 
polarizability (3 of the proton so small (/3 P ~ 2 x 10~ 4 fm 3 despite of a very large 



paramagnetic contribution of the A resonance which, in various evaluations based on quark 
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models, dispersion theories, effective Lagrangians, etc., ranges between +7 to +13 x 10 4 fm 3 
(see, for example, j55,56|). In this context, it is sometimes stated that the pion cloud pro- 



duces a large diamagnetic (i.e. negative) contribution to (3 owing to the Langevin mechanism 
producing a negative magnetic susceptibility of bound charges. This point of view became 
especially popular after calculations of Pn in the Skyrme model (see, for example p7|-p9[|), 
in which the pion field of the soliton produces indeed a very big diamagnetic susceptibil- 
ity resulting from the pion-cloud tail of the soliton, ranging from —8 to —16 x 10~ 4 fm 3 
|57^p9[. Figure [5], however, shows that the pion-cloud periphery with r > 1 fm carries a 
very small diamagnetism of only —0.45 x 10~ 4 fm 3 , where the last number is expected to be 
a reliable estimate because ChPT, even at leading nontrivial order, should be reasonable at 
such distances. We conclude that Skyrme models overestimate the pion-cloud contribution 
to diamagnetism and that a source for the missing diamagnetism is probably related with 
short-range mechanisms rather than with the pion cloud itself. In some dispersion theories 
of Compton scattering ||60|| , an additional exchange with a hypothetical cr-meson is invoked 
in order to generate agreement with existing experimental data. This, of course, is just an 
oversimplified model for such a short-range contribution. 

It would be very interesting and instructive to extend the presently available ChPT pre- 
dictions for the generalized polarizabilities of the nucleon at least to order 0{p i ). This would 
allow one to check whether a modification of the density /3(r), including relativistic (nucleon 
recoil) effects and other higher-order corrections, indeed provides sufficient diamagnetism 



VI. SUMMARY AND CONCLUSIONS 



In the present paper we have developed a covariant formalism leading to a parametriza- 
tion of the VCS tensor in terms of Lorentz invariants free from kinematical singularities and 
constraints. We motivated and performed a gauge-invariant division of the VCS tensor into 
contributions of generalized Born terms and a structure-dependent residual part. We then 
discussed the case of a real final photon and a space-like virtual initial photon which can be 
described in terms of three invariant amplitudes depending on three kinematical variables. 
In the limit q J — >• 0, the three residual amplitudes reduce to functions of q 2 only which we 
identified as generalized dipole polarizabilities a^g 2 ), ar(q 2 ), and f3(q 2 ). All of them can, in 
principle, be determined in virtual Compton scattering, although the transverse electric po- 
larizability is not accessible in experiments sensitive to structure-dependent effects of 0(q') 
only. We proposed a physical interpretation of these polarizabilities in terms of spatial dis- 
tributions of an induced electric polarization and magnetization, respectively. In particular, 
we argued that a knowledge of all three polarizabilities is required for a full description of 
induced polarization phenomena. Following this line, we calculated spatial distributions for 
pions, kaons, and the baryon octet as Fourier integrals, using the predictions of ChPT. It 
was found that the distributions obtained confirmed expectations based on a picture of a 
hadron's periphery caused by a "classical" pion cloud. Of course, any practical analysis 
of experimental data on photon scattering will eventually deal with the original, precisely 
defined momentum-space form factors, polarizabilities, etc. Thus, one should handle the 
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found spatial distributions with care. On the other hand, the r-space interpretation of such 
g-dependent quantities clearly allows for a more intuitive visualization in analogy to the 
phenomenology and terminology of classical electrodynamics and nonrelativistic quantum 
mechanics. 
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APPENDIX A: SPATIAL DISTRIBUTIONS TO ONE LOOP 



In this appendix we collect the analytical results for the spatial distributions in ChPT 
at the one- loop level which are easily obtained through Eq. (|5.4j ). In the formulas below we 
use the notation 

t = q 2 , x w = — , x K = —3-, z w = 2m 7r r, z K = 2m K r. (Al) 



m; m 



K 



The imaginary parts of the scalar and vector form factors of the pion, and of the isovec- 
tor charge form factor of the nucleon are determined by the function J^\x) which has a 
branching point at x = 4: 

2r — 1 

ImF v (t) = ^Av(x w ), 
o 



where 



ImG£(t)=(^— + )V(x,), (A2) 



Evaluating the integral representation of Eq. (|5.4| ) with the above imaginary parts results in 
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Anr 2 F s (r) 



Anr 2 F v (r) 



4^r 2 G v E (r) 



mlF s (0) 
(4?rF w ) 2 


r 48 r ^ , 




g 

— 






8 + 40^1 


~ (4tt^) 2 I 


^7T 



96 + 14z 2 



16 



(A4) 



where i£z,(z) is the modified Bessel function, K v (z) = / °° e zcosht cosh(z4) (it. 



The generalized polarizabilities of the nucleon, Eq. (|5.14|) , at time-like momenta q have 
both a cut starting at t — 4m 2 and a pole at i = 4m\: 



7T 



ImQi(t) = — E. 



Am 2 \ 6 Ht - Am 2 ) „, 9N 

2 + ^ ) + 4m, 5(t - 4ml) 



7T 



Im/3(f) = — El 



t > Vi 

4m 2 \ 6>(t -4m 2 ) 



t 



- 4m n 5(t - 4m 2 n ) 



(A5) 



Accordingly, Eq. (|5.4j) results in 



El 



4nr 2 a L (r) = [(1 + z^eT** + -f Ei(-z n ) 
4nr 2 f3(r) = 3- [(1 - z^e"** - ^ Ei(-^) 



(A6) 



where Ei(— z) = — J^°{e t /t) dt is the exponential integral. The contribution of kaon loops 
is obtained by the substitutions E n — > Ek and m n — > mx as explained in Section M. 



For pions and kaons, the generalized polarizabilities of Eq. ( |5. 19| ) have branching points 
at t — 4m 2 and 4m 2 K but no poles: 



lm(3(t) 



2tt 



4ttM (47TF,) 2 
Then Eq. (O) gives 



x * <lxl- 4x„ 



x 2 K — 4xk 



(A7) 



4irr 2 /3(r) 



2m n C 7T z 7T K Q (z n ) + 2m K C K ZKK (z K ) 



4nM (4ttF,) 2 

+ -^B n zl(K 1 (z 7T )- J z K (x)dx 



(A8) 
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